Continuous variable entanglement of phase locked light beams 
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We explore in detail the possibility of intracavity generation of continuous- variable (CV) entangled 
states of light beams under mode phase-locked conditions. We show that such quantum states can 
be generated in self-phase locked nondegenerate optical parametric oscillator (NOPO) based on a 
type-II phase-matched down-conversion combined with linear mixer of two orthogonally polarized 
modes of the subharmonics in a cavity. A quantum theory of this device, recently realized in the 
experiment, is developed for both sub-threshold and above-threshold operational regimes. We show 
that the system providing high level phase coherence between two generated modes, unlike to the 
ordinary NOPO, also exhibits different types of quantum correlations between photon numbers and 
phases of these modes. We quantify the CV entanglement as two-mode squeezing and show that 
the maximal degree of the integral two-mode squeezing(that is 50% relative to the level of vacuum 
fluctuations) is achieved at the pump field intensity close to the generation threshold of self-phase 
locked NOPO, provided that the constant of linear coupling between the two polarizations is much 
less than the mode detunings. The peculiarities of CV entanglement for the case of unitary, non- 
dissipative dynamics of the system under consideration is also cleared up. 

PACS numbers: 03.67.Mn, 42.50.Dv 
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I. INTRODUCTION 

It is now believed that entanglement of quantum com- 
posite systems with continuous degree of freedom is the 
basis of most apphcations in the field of Quantum In- 
formation [l|. Interest in continuous- variable (CV) en- 
tanglement is being extensively excited by successful ex- 
periments on quantum teleportation based on two-mode 
squeezed states i2| as well as the experiments dealing with 
entanglement in atomic ensembles j^l- Since then, re- 
markable theoretical and experimental efforts have been 
devoted to generating and quantifying CV entangled 
states. 

In this paper we propose a novel type of CV entan- 
gled states of light-field with reduced phase noise. They 
are different from the well-known entangled Einstein- 
Podolsky-Rosen (EPR) states generated in a nondegen- 
erate optical amplifier 4, 5], which exhibit large phase 
fluctuations. We believe that such entangled states of 
light-field with localized phases can be generated in a self- 
phase locked nondegenerate optical parametric oscilla- 
tor (NOPO), based on the type-II phase-matched down- 
conversion and additional phase locking process stipu- 
lated by the intracavity waveplate. The motivations for 
this study are the following: 

For the first time the CV entangled states of light 
were studied in [J| and demonstrated experimentally 
in Q for nondegenerate optical parametric amplifier 
(NOPA). Then a CV entanglement source was built from 
two single-mode squeezed vacuum states combined on 
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a beamsplitter l3|. It is well known that each of the 
orthogonally polarized and frequency degenerate fields 
generated by NOPO is a field of zero-mean values. The 
phase sum of generated modes is fixed by the phase of 
the pump laser, while their phase difference undergoes 
a phase diffusion process ^ stipulated by vacuum fluc- 
tuations. As a rule, the NOPO phase diffusion noise 
is substantially greater than the shot noise level, that 
limits the usage of NOPO in precision phase-sensitive 
measurements. Various methods based on phase locking 
mechanisms 0, B H UOl have been proposed for re- 
ducing such phase diffusion. In the comparatively simple 
scheme realized in the experiment |6|, self- phase locking 
was achieved in NOPO by adding an intracavity quarter- 
wave plate to provide polarization mixing between two 
orthogonally polarized modes of the subharmonics. The 
evidence of self-phase locking was provided there by the 
high level of phase coherence between the signal and idler 
fields. Following this experiment, the semiclassical the- 
ory of such NOPO was developed in Q- Recently, the 
schemes of multiphoton parametric oscillators based on 
cascaded down-conversion processes in x media placed 
inside the same cavitv and showing self-phase locking 
have been proposed [g. As was demonstrated in Q, the 
system based on combination of OPO and second har- 
monic generation also displays self-phase locking. The 
formation of self-phase locking and its connection with 
squeezing in the parametric four-wave mixing under two 
laser fields has been demonstrated in [lOJ. An impor- 
tant characteristic of self-phase locked devices concerns 
the phase structure of generated subharmonics. Indeed, 
the formation of the variety of distinct phase states un- 
der self-phase locked conditions has been obtained in the 
mentioned Refs. 0, S 0; U^- It was recently noted 
that the schemes involving phase locking are potentially 



useful for precise interferometric measurements and op- 
tical frequency division because they combine fine tun- 
ing capability and stability of type - II phase matching 
with effective suppression of phase noise. That is why 
we believe it will be interesting to consider phase locked 
dynamics also from the perspective of Quantum Optics 
and, in particular, from the standpoint of production of 
CV entanglement. 

A further motivation for such task is connected with 
the problem of experimental generation of bright entan- 
gled light. So far, to the best of our knowledge, there 
is no experimental demonstration of CV entanglement 
above the threshold of NOPO. The progress in experi- 
mental study of bright two-mode entangled state from 
cw nondegenerate optical parametric amplifier has been 
made in [Ij • The theoretical investigation of CV entan- 
gled light in transition through the generation threshold 
of NOPO is given in [l3| . One of the principal experimen- 
tal difficulties in advance toward a high-intensity level is 
the impossibility to control the frequency degeneration 
of modes above the threshold. We hope that the usage 
of phase locked NOPO may open a new interesting pos- 
sibility to avoid this difficulty. 

In this paper we report what is believed to be the first 
investigation of self- phase locked CV entangled states. 
We develop the quantum theory of self-phase locked 
NOPO, with decoherence included, in application to the 
generation of such entangled states. This scheme is based 
on the combination of two processes, namely, type-II 
parametric down-conversion and linear polarization mix- 
ing with cavity-induced feedback. The parametric down- 
conversion is a standard technique used to produce an 
entangled photon pairs as well as CV two-mode squeezed 
states U • The beam splitter including polarization mixer 
are also considered as experimentally accessible devices 
for production of entangled Ught-fields [131 . Besides 
these, there have been some studies of a beam splitter 
for various nonclassical input states, including two-mode 
squeezing states [ij. It is obvious, and also follows from 
the results of the mentioned papers [^ |^ [13 , that the 
operational regimes of the combined system with cavity- 
induced feedback and dissipation drastically differ from 
those for pure processes. We show below that analogous 
situation takes place in the investigation of quantum- 
statistical properties of a combined system such as the 
self-phase locked NOPO. 

The paper is arranged as follows. In Section II we 
formulate the model of combined NOPO based on the 
processes of two-photon splitting and polarization mix- 
ing, and present a semiclassical analysis of the system. 
Section III is devoted to the analysis of quantum fluctu- 
ations of both modes within the framework of lineariza- 
tion procedure around the stable steady-state. In Section 
IV we investigate the CV entangling resources of self- 
phase locked NOPO on the base of two-mode squeez- 
ing for both sub-threshold and above-threshold opera- 
tional regimes. We also discuss in Section IV the case 
of unitary, non-dissipative dynamics using a well justi- 



fied small-interaction time approximation. We summa- 
rize our results in Section V. 



II. MODEL OF SELF-PHASE LOCKED NOPO 

As an entangler we consider the combination of two 
processes in a triply resonant cavity, namely, the type - 
II parametric down-conversion in x'^''- medium and po- 
larization mixing between subharmonics in lossless sym- 
metric quarter-wave plate. The Hamiltonian describing 
intracavity interactions is 

+ihk (e^'^^ay,a^a2 — e^'^'^''at^aia2) 
+ftx(e'*^afa2 + e-**-aia+), (1) 

where a^ are the boson operators for the cavity modes 
LOi. The mode 03 at frequency uj is driven by an exter- 
nal field with amplitude E and phase $l, while ai and 
02 describe subharmonics of two orthogonal polarizations 
at degenerate frequencies lli/2. The constant A:e'*'= deter- 
mines the efficiency of the down-conversion process. Lin- 
ear coupling constant denoted as xe**'' describes the en- 
ergy exchange between only the subharmonic modes and 
besides, x is determined by the amount of polarization 
rotation due to the intracavity waveplate, <^j^ determines 
the phase difference between the transformed modes. We 
take into account the detunings of subharmonics A^ and 
the cavity damping rates 7^ and consider the case of high 
cavity losses for pump mode (73 ^ 71,72), when this 
mode can be adiabatically eliminated. However, in our 
analysis we take into account the pump depletion effects. 
The principal scheme of this device, the so-called self- 
phase locked NOPO, is shown on Fig. ^ for the special 
configuration when the coupling of in- and out-radiation 
fields occurs at one of the ring cavity mirrors. Note, that 
the comparison of this configuration with the analogous 
one considered in |7| , where the pump field is a travelling 
wave, is discussed below. 

Due to explicit presence of dissipation in this problem, 
one has to write the master equation for the reduced den- 
sity matrix of the system. The reduced density operator 
p within the framework of the rotating wave approxima- 
tion and in the interaction picture is governed by the 
master equation 



dp 
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'a^) 



(2) 



where 



H' = ^ hAiufai + ihE (aj - 03) 

4=1 

-I- ihk{a3,al a2 — a;^aia2) + hx{aia2 + aia^).(3) 




FIG. 1: The principal scheme of self-phase locked NOPO in 
triply resonant cavity for the pump mode at frequency u and 
two subharmonic modes of two orthogonal polarizations {]) 
and (^) at frequency uj/2. We have chosen the special path 
for the pump mode to underline that the pump mode rapidly 
decays (73 3> 71,72) and is eliminated adiabatically. The 
type-II phase-matching condition for the process oj ^ ^ (t) + 
^ (^) is satisfied in x medium. The intracavity quarter 
waveplate QWP provides a polarization mixing between the 
orthogonally polarized subharmonics. 



Let us also note that this equation is rewritten through 
the transformed boson operators ai —^ aiexp(— i$i), 
with $, being $3 = $l, $2 = 5 ($l -H $fe - $x)' 
•J*! = 5 i^L + ^k + $x)- That leads to cancehation of 
the phases on the intermediate stages of calculations. As 
a result, the Hamiltonian H' depends only on real-valued 
coupling constants. 

In order to proceed further, we now consider the 
phase-space symmetry properties of the two subharmonic 
modes. It is easy to check that the interaction Hamil- 
tonian satisfies the commutation relation [H' , U (tt)] = 



with the operator U (9) — exp [i6 (a^ai — 0^02)] . More- 
over, analogous symmetry [p (i) , U (tt)] ~ takes place 
for the density operator of the system, which obeys the 
master equation (0). Using this symmetry we establish 
the following selection rules for the normal-ordered mo- 
ments of the mode operators of the subharmonics: 



{a{; aj a^a^) = 0, 



(4) 



iik + l + m + n^ 2j, j = 0, ±1, ±2, .... Another pe- 
culiarity of NOPO with Hamiltonian ^ is displayed in 
the phase space of each of the subharmonic modes. The 
reduced density operator for each of the modes is con- 
structed from the density operator p by tracing over the 
other modes pi = Tr2Tr3 (p), p2 = TriTr^ (p). There- 
fore, we find that [pi (t) , C/i (tt)] = [p2 (t) , U2 (tt)] = 0, 
where Ui (9) = exp(^i9afai) , (i = 1,2) is the rotation 
operator. It is easy to check using these equations that 
the Wigner functions Wi and W2 of the modes have a 
two-fold symmetry under the rotation of the phase-space 
by angle tt around its origin. 



W^ (r, 9) = W, (r, 9 + tt) , 



(5) 



where r, 9 are the polar coordinates of the complex phase 
space. Note, that such symmetry relationships Q, (O 
radically differ from those taking place for the usual 



NOPO without the quarter-wave plate. Indeed, in this 
case, (case of x = in Eqs. CQl-®) the Wigner functions 
Wi are rotationally symmetric and the symmetry prop- 



erties Q read as (a^ 



+kA „+m 



a\a 



1"2 



a^) = 0, if fc - Z 9^ 



We perform concrete calculations in the positive P- 
representation 15] in the frame of stochastic equations 
for the complex c-number variables a^ and Pi correspond- 
ing to the operators ai and af 

-K7- = - i'Ji + i^i) ai + ka3l32 - ix<^2 + Ri, (6) 
^ = ~i-f,-iAi)Pi + kP3a2+ixP2 + Rt, (7) 



-7— ^ --/sas + E - kaia2, 
ot 



(8) 
(9) 



The equations for a2 , (32 are obtained from lO ) , O 
by exchanging the subscripts (1)«=^(2). R12 are Gaussian 
noise terms with zero means and the following nonzero 
correlators: 



{Ri{t)R2{t')) ^ ka^S {t - t') , 



{R+{t)R+{t'))^kM{t-t') 



(10) 



(11) 



Recall, that we consider the regime of adiabatic elimina- 
tion of the pump mode. In this approach the stochastic 
amplitude 03 and 13^ are given by 



az = [E ~ kaia2) hz- 



p3 ^ (E ^ kl3il32) /j3- 



(12) 



(13) 



Substituting amplitudes a^, (3^ into Eqs.®, and l(Tn|l . 
(|ll|l we arrive at 



•9"! / AX 

^^ = - (71 +iAi)ai 



+ (e - \a1a2) (32 - ixct2 + Ri 



(14) 



^ = -(7i-zAi)/3i 



+ {e-Xpi(32)a2+ixP2 + Rt- (15) 



{Ri{t)R2{t')) = (e - Xaia2) 6 {t - t') , (16) 



{R+{t)R+{t')) - (e - A/3i/32) 5 (t - t') . (17) 

Here £ = kE/js, A = fc^/73. So the Eqs.lUH), IjTsl) and 
corresponding Eqs. for a2 , (32 involve the depletion effect 
of the pump mode, which leads to the appearance of the 
above-threshold operational regime. 



It should be noted that Eqs.JSl, ^ for pump mode 
ampHtude do not involve terms of the linear coupling 
between the modes of subharmonics. From this we can 
immediately recognize that the Eas. (|12|l . Ijl^l) and hence 
Eas. H14|l . (|15|1 take place for arbitrary values of the pa- 
rameter X- Nevertheless, the effect of linear coupling be- 
tween the modes are displayed in the dynamics of pump 
mode through the amplitudes ai and a2- In the adi- 
abatic approach these effects are described through the 
pump depletion terms kaia2/j3, fc/3i/32/73 in the expres- 
sions for the pump field amplitudes H12|l . H13|) and hence 
are completely taken into consideration in Eas. ll4|l . 115|) 
as well as in the correlators ((TB)l . ((TTJI . 

First, we shall study the steady-state solution of the 
stochastic equations in semiclassical treatment, ignoring 
the noise terms for the mean photon numbers rijo and 
phases ipjo of the modes {uj = ajf3j, (pj = 57 In {aj/(3j)). 
The mean photon numbers read as 



nin = ^ 



A2 
Ai 






1/2 r 



1/2 



(£?.)' +7^ 



+ 7—7 



where 



7 = 



71 f ^2 

Ai 



1/2 



72 / Ai 
2 U2 



1/2 



(18) 



(19) 



(20) 



ordinary NOPO. These results are in accordance with 
the ones obtained in [y, Q , but for another configuration 
of NOPO. In the scheme proposed in [g only the signal 
and idler modes are excited in the cavity, while the pump 
field is a travelling wave. Nevertheless, in the adiabatic 
regime considered here there is correspondence between 
both schemes. Indeed, it is not difficult to check that the 
results (|18f) . H19|) transform to the corresponding results 
of the mentioned scheme [3 by replacing the parameter 
£ with the corresponding pump field amplitude. 

We now turn to the standard linear stability analysis of 
these solutions, assuming for simplicity, the perfect sym- 
metry between the modes, provided that the cavity decay 
rates and the detunings do not depend on the polariza- 
tion ( 71 = 72 = 7, Ai = A2 = A). The stability of the 
system is governed by the matrices F and F^ , F^ describ- 
ing the dynamics of small deviations Sai and 6/3i from 
the semiclassical steady-state solutions (see. Sec. III). We 
reach stability if the real parts of eigenvalues of these 
matrices are positive. This analysis displays an evident 
dependence on the sign of the detunings Ai = A2 = A. 
For the positive detuning, A > 0, only the steady-state 
solutions nj~Q = n^g and (/s^q — ip^Q, (/s^q -|- ip^^ are stable, 
while for the case of negative detuning the stability holds 
for the solutions with the (-I-) superscript. As this anal- 
ysis shows, for either sign of detuning, the threshold is 
reached at e ^ eth, where 



£th = V (x - U 



+ r, 



(25) 



As we can see from (jl8|l . H19|l . nj^ and njo are real and 
positive for e exceeding two critical points 

{etf = 7172 + A1A2+X' 



and the steady-state stable solution for mean photon 
numbers can be written in the general form as 



1 

no = nw = n2o = - 



e'^-{x-\A\f-j 



(26) 



TV4x'AlA2-(7lA2-72Al)^ (21) The phases are found to be 



and, besides, the solutions nj, and n,n, (i=l,2) corre- 



spond to two distinct values e+, and e^^ accordingly. The 
steady-state values of the phases corresponding to each 
of the critical points e^^, e^ are obtained as 



sin(99. 



20 



fw) = sin(^20 - fw) 



1 

2^ 



71 



A2 
Ai 



1/2 



72 



Ai 
A2 



1/2' 



(22) 



cos(^± + pfo) - iv'e2_(,±)%^2. (23) 

It is easy to check that these solutions exist for both 
modes only if the following relation holds 



4x'AiA2> (71A2-72A1) 



(24) 



Let us note that the steady-state solutions (|22|l , H23|) com- 
pletely determine the absolute phases of the orthogonally 
polarized modes, which are hence self-locked, unlike the 



<^io = y'20 = --Arcsin - (x + |A|) + nk, (27) 

2 £ 

for A > 0. For the opposite sign of the detuning, A < 0, 
as we noted, the mean photon numbers are given by the 
same Ea. (l26|l . while the phases read as 

<Pio = ^Arcsm-(x+\A\)+Tr(k + -y (28) 
<P20 = -Arcsin-(x+|A|)+7r(^fc--j, 

(A: = 0, 1, 2, ..). In the region e ^ eth the stability condi- 
tion is fulfilled only for the zero amplitude steady-state 
solution ai ^ a2 ~ (3i — f32 — 0. So, the set of above- 
threshold stable solutions for both modes have two-fold 
symmetry in the phase-space which was indeed expected 
from symmetry arguments ©. 

Let us now consider the output behavior of the sys- 
tem for the special scheme of generation, when the cou- 
plings of in- and out-fields occur at only one of the ring- 
cavity mirrors (see Fig. ^. Taking into account that 



only the fundamental mode is coherently driven by the 
external field with (a™) = ^/73, while the subharmonic 
modes are initially in the vacuum state, we obtain for 
the mean photon numbers (in units of photon number 
per unit time) n|" = E'^/2-fi, ?i°"* = 27^7^,0 (i = 1,2) 
and hence 71°"* — 7^2"* — 27710. Accordingly, parametric 
oscillation can occur above the threshold pump power 
Pth = fiiu^E^i^/2j3, where the threshold value of the 
pump field is equal to 



where /i — 1,2,3,4 and Sa^ = {Sai,Sa2,Sa3,Sa4) = 
{Sai,5a2,6l3i,Sf32), R^ = {Ri,R2,Rt,Rt)- The 4x4 
matrix Ffj,i, is written in the block form 



F 



with 2x2 matrices 



A, B 
B\ A* 



(31) 



i?*/. = fV(x-|A|)'+72. 



(29) 



We are now in a position to study quantum effects in 
self-phase locked NOPO and will state the main results 
of the paper concerning CV entanglement. 



III. ANALYSIS OF QUANTUM 
FLUCTUATIONS 



«X, 7 + 7^/' V 1' 



The noise correlators are determined as 



(32) 



(i?^ (a, t) i?" (a, t')) = D^, {a) 5 {t ~ t') (33) 
with the following diffusion matrix 



The aim of the present section is to study the quantum- 
statistical properties of self-phase locked NOPO in linear 
treatment of quantum fluctuations. Quantum analysis 
of the system using P- representation is standard. A 
detailed description of the method can be found in [T^ ■ 
We assume that the quantum fluctuations are sufficiently 
small so that Eas. (|14|l . (|15|l can be linearized around the 
stable semiclassical steady state ai (t) = a° -I- 5ai (i), 
Pi (i) — Pi + SPi (t) . This is appropriate for analyz- 
ing the quantum-statistical effects, namely CV entan- 
glement, for all operational regimes with the exception 
of the vicinity of threshold, where the level of quantum 
noise increases substantially. It should also be empha- 
sized at this stage that in the above-threshold regime of 
self-phase locked NOPO the steady-state phases of each 
of the modes are well-defined in contrast to what hap- 
pens in the case of ordinary NOPO, where phase diffusion 
takes place. According to this effect, the difference be- 
tween the phases, as well as each of the phases, can not 
be defined in the above-threshold regime of generation 
of the ordinary NOPO. On the whole, the well founded 
linearization procedure can not be applied for this case. 
Nevertheless, the linearization procedure and analysis of 
quantum fluctuations for ordinary NOPO become pos- 
sible due to the additional assumptions about temporal 
behavior of the difference between the phases of the gen- 
erated modes jTq . 

We begin with consideration of below-threshold oper- 
ational regime, for E < Eth, where the equations lin- 
earized around the zero-amplitude solution can be writ- 
ten in the following matrix form 



^^6a^ ^ -F^Ja" + R^ ia,t) , 



D = 



B, 

0, B* 



(34) 



First, we calculate the temporal correlation functions 
of the fluctuation operators. The expectation values of 
interest can be written as the integral 



(mO 



fiiy 



/niinir,r i 
dT (e^(*--)De^"(*'-^)) 
-00 

(35) 
F"^ being the transposition of the matrix F. The integra- 
tion over dr can be performed using the following useful 
formula for operators DF^ = FD, obtained by straight- 
forward calculation. As a consequence, we arrive at the 
expression De^ * — e^*'D. Finally we obtain 



1 



{5a^[t)5a'^{t')) = --[F e 



and hence for t = t' 

{Sa^" {t) ia" it)) 



i^-f|*-*'I 



D 



(36) 



-{F-^D) , 



(37) 



This formula, however, is not very convenient for prac- 
tical calculations. Therefore, we rewrite the correlation 
functions of the quantum fluctuations ^a^, bPi in a more 
simple form through the two-dimensional column vectors 
(5a — ((5ai,(5a2) , <5/3 = ((5/3i,(5/32) ■ Upon performing 



(30) the calculation, we finally arrive at 



da {6a) ) ~ 



5a{5py 



£ 



5-4 _ 4A2;^2 
-2 



£ 



2(S'4-4A2x2) 



-2xA, 52 X 
S\ -2xA j 

^', -2xA 

-2xA, 52 



r 



X(52-2A2), A (52 -2x2 
A (52 _ 2x2), x(52-2A2 



(38) 



where 5'^ is introduced as 



52 



7^ + x' + A 



2 £2_ 



(39) 



Note, that S'^ > in the below-threshold regime. 

As an apphcation of these results the mean photon 
number in the below-threshold regime can be calculated 
as follows. 



2 c2 



m = 712 



e'S- 



2 (S'4 - 4A2x2) ■ 



(40) 



Next we focus on the mode locked regime, for E > Eth, 
considering the Eqs. p4|) . (|15|l in terms of the fluctua- 
tions Srii (t) = Hi (t) — riio and dipt (i) = (pi (t) — tpio 
of photon number and phase variables. In this regime 
the dynamics described by the linearized equations of 
motion actually decouples into two independent dynam- 
ics for two groups of combinations 6n± = Sn2 ± 6ni, 
5ip± — Sip2 ± 6ipi. In fact, one has 



dt \6ip+, 



dt\Sip- 



= -F_ 



6n^ 
Sip^ 

Sri- 

Sip- 



Rn-\ 



Rn- 
Rio- 



(41) 



(42) 



The drift, F±, and the diffusion matrices, 
{R.{t)R,{t')) = D^+^6it-t'), (z,j = 7i+,5^+); 

{Rn{t)Rrn{t')) = D^S {t - t') , {n, m = n^,ip-) are 
respectively 



F.= 



2Xno, 4noesin((p20 + '/'lo) 
0, 2 (7 + Ano) 



(43) 



^ / 27, 4noXsinn/A 
^ A/no, 



(44) 



I?(±) = ± 



4no7, - 2e sin(^20 + y^w) 
-2£ sin(v320 + V'lo), - ^/no 



(45) 



Due to the decoupling between (+) and (— ) 
combinations of the modes, we conclude that the 
following temporal correlation functions are equal 
to zero, {Sip+ (t) Sip^ {t')) = {Sti+ (t) 6n_ {t')) = 

{dn± {t)dip^ (t')) = 0, and hence ({Sipi) ) = ({6ip2) ), 

({6ni) \ = ((6712) ). The other correlation fmictions 

can be calculated in the same way as described above for 
the below-threshold regime. The temporal correlation 
functions are derived as 



G;; Si <*"*<''> '*>*<'■» 



f-ie-'i'l'-'ICi, 



and hence 



6n^ 
S(pj 



{Sn±,S^±)) = --F^'D, 



(46) 



(47) 



Performing the concrete calculations for each of the 
cases A > and A < 0, we arrive at the following results 



Sip+ 



{Sn+,6ip+) 
((5ri_, dip-) 



4^0 [7(7 + Ano) + (x- I A 1)2] 
4Ano(7 + Ano) V -2Ano(x- | A \)sign{A), 



1 



-Ano(x- I A \)sign{A) 
-A7 



4riox(x- I A I), 2x"fsign{A) 

4|A|xV 2x7si5"(A), ^{j^- \ A \ {x~ \ A\) 



(48) 
(49) 



r 



We see that the considered system displays different 
types of quantum correlations in terms of the stochas- 
tic variables, namely, between photon- number sum and 



phase sum in the modes, as well as between photon- 
number difference and phase difference in the modes. 
These results radically differ from the ones taking place 



for usual NOPO, where the correlation between photon- 
number difference and phase sum in the modes is realized. 
The results obtained indicate the possibilities to produce 
entanglement with respect to the new types of quantum 
correlations. 



IV. CV ENTANGLEMENT IN THE PRESENCE 
OF PHASE LOCKING 



Let us now turn our attention to quantum statisti- 
cal effects and the entanglement production for the case 
of perfect symmetry between the modes (71 =72 =7, 
Ai = A2 = A). Our aim in this section is to study the 
interplay between phase locking phenomena and CV en- 
tanglement for the self-phase locked NOPO. We note that 
unlike the two-mode squeezed vacuum state, the state 
generated in the above-threshold regime of NOPO is non- 
Gaussian, i.e. its Wigncr function is non-Gaussian |l7l |. 
Recently, it has been demonsrated [Tg, that some sys- 
tems involving beam splitters also generate non-Gaussian 
states. The general consideration of this problem for 
self-phase locked NOPO seems to be very complicated. 
However, the mentioned results allow us to conclude 
that the state generated in self-phase locked NOPO is 
most probably non-Gaussian. So far, the inseparabil- 
ity problem for bipartite non-Gaussian state is far from 
being understood. On the theoretical side, the neces- 
sary and sufficient conditions for the separability of bi- 
partite CV systems have been fully developed only for 
Gaussian states, which are completely characterized by 
their first and second moments. To characterize the 
CV entanglement we address to both the inseparabil- 
ity and strong EPR entanglement criteria [13 which 
could be quantified by analyzing the variances of the 
relevant distance V- = V {Xi — X2) and the total mo- 
mentum Vj^ = V [Yi + Y2) of the quadrature amplitudes 
of two modes Xk — -^ [a^ exp (— 16'^) -I- a^ exp (i^^)] , 

>fc = ^ [at cxp (-iOk) - Ofe exp {i6k)\ , (fc = 1, 2), where 

V{X) = {X'^) — (X) is a denotation for the variance and 
61, is the phase of local oscillator for the k-th mode. The 
two quadratures Xk and Yk are non commuting observ- 
ables. The inseparability criterion for the quantum state 
of two optical modes reads as jl^ 



V 



'-{V+ + V-) < 1, 



(50) 



i.e. indicates that the sum of variances drops below the 
level of vacuum fluctuations. Since the states of the sys- 
tem considered are non-Gaussian, the criterion (|50|l is 
only sufficient for inseparability. The strong EPR entan- 
glement criterion is quantified by the product of variances 
as V+y_ < J. We remind that the sufficient condition for 
inseparability 15U|1 in terms of the product of variances 
reads as V+y_ < 1, i.e. is weaker than the strong EPR 
condition. 

In order to obtain the general expressions for variances, 
we first write them in terms of the boson operators corre- 



sponding to the Hamiltonian QJ . We perform the trans- 
formations ai — > Qi exp (z<&i), which restore the previous 
phase structure of the intracavity interaction. Using also 
the symmetry relationships (@J) we find quite generally 
the variances at some arbitrary quadrature phase angles 
01, 02 as 



where 



V± ^V±Rcos{A0), 

-2 I (aiaa) | cos(I]6l + $arg), 
ii; = 2Re((a?)e*^'')-2| (a+a^) 



(51) 



(52) 



(53) 



A0 = 02 - 01 - $x> S0 = 01 + 02 + *( + $fc, (54) 

and n = (oj^ai) = (0^02) is the mean photon num- 
ber of the modes, $arg — arg(aia2). So, in accordance 
with the formula (|51|l , the relative phase <i>^ between the 
transformed modes gives the effect of the rotation of the 
quadrature amplitudes angle 02 — ^i ■ 

Obviously, the variances V± and hence the level of CV 
entanglement depend on all parameters of the system 
including the phases. The minimal possible level of V 
is realized for an appropriate selection of the phases 0^, 
namely for 0i -I- 02 = — arg (0102) — $; — $fc. Further, in 
most cases we assume that this phase relationship takes 
place, but do not introduce new denotations for V and V± 
for the sake of simplicity. In this case, (S0 = — 'I'arg), in 
correspondence with the formula (|51(l . the variances V± 
depend only on the difference between phases A0 and we 
arrive at 

V = l + 2{n-\{aia2)\), (55) 

i? = 2 Rc((a?) e-**"s ) - 2 | (0+02) | . (56) 

For the NOPO without additional polarization mixing 
(af a2) ~ (af) = (a^) = and hence the case of the 
symmetric variances V+ = V- is realized. The phase- 
locked NOPO generally has non symmetric uncertainty 
region. However, the variances V+ and V- become equal 
for the special case of 02 — 0i — <&x ^ f ' when the insep- 
arability condition reads as V < 1, {V = V- =V+). We 
note, that the relative phase <I>^ plays an important role 
in specification of the entanglement. 



Entanglement in the case of unitary 
time-evolution 



So far, we have considered mainly steady-state regime 
of generation, including the effects of dissipation and 
cavity-induced feedback. However, in order to better 
understand the pecularities of the entanglement for the 



system under consideration, it would be interesting and 
desirable to study the situation, when the dissipation in 
the cavity is unessential and the evolution of the modes, 
due to their interactions, is described approximately by 
the master equation jSJ without losses. In our analysis 
we shall assume that both subharmonic modes are ini- 
tially in the vacuum state. Then, the system's state in 
the absence of losses is generated from the vacuum state 
by the unitary transformation 



|*(i))=u(t)|0)j0)2-e 



_ ^-fHi-att 



o)i 10)2 



(57) 



where 



Hi, 



ih- 



kE 
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aia2) + hx {a^a2 - fliaj) . (58) 



Here we neglect the pump depletion effects considering 
classically the amplitude of driving field as a constant 
and also consider the case of zero-detunings Ai = A2 = 
A = 0. For the cavity configuration considered the va- 
lidity of such approximation is guaranteed for short in- 
teraction time intervals {kE/^3)~^ ,x~^ £ ^* ^ 7~^7 
provided that the coupling constants kEf-f^ and x ex- 
ceed the dumping rates of the modes. 

Note, that for x — > the state (|57ll is transformed 
to the well-known two mode squeezed state 20]. We 
note, that while the entangled two-mode squeezed vac- 
uum state is achieved in the laboratory (in the nonde- 
generate parametric amplifier [2|), the generation of the 
state (|57|l is an interesting option for the future. From an 
experimental point of view it can be realized in the com- 
bined scheme, that is parametric down conversion and 
wave plate, driven by a pulsed pump field. It is obvious, 
that the interaction time Ai in this case is limited by the 
duration T of laser pulses, Ai < T. Roughly speaking, 
such scheme is similar to the optical nondegenerate para- 
metric amplifier in the presence of phase-locking process, 
where the interaction time is increased by a cavity for 
both subharmonic modes. However, we note again, that 
a complete study of the experimental configuration re- 
lated to the nondissipative case is beyond the scope of 
this paper. 

Now let us examine the entangled properties of the 
resultant state H57|l. Using formula (I52II . where n — 
{u^^ {t)af aiu{t)) and (ai(i)a2(t)) = {u^^ {t)aia2u{t)) , 
after a long but straightforward algebra we get the final 
result for the variance V = ^(Vl). + V-) for two different 
operational regimes: £ < x. ^^id e > x- The variance for 
the range e < x of 3' comparatively weak scaled pump 
field e = kEj^-^ reads as 



V[i) 



[cos{2fit) - 1] 



- sin(2/xi) cos (116') , 



(59) 



where /i — \/x^ — e^. Thus, we observe a periodic evo- 
lution of V in this regime typical for the linear coupling. 




FIG. 2: Unitary evolution of the variance V (t) versus the 
scaled dimensionless interaction time et for the range of e < x 
and provided that cos {T,9) — 1. The parameters are: e/x = 
0.2 (curve 1), e/x ~ 0.4 (curve 2) and e/x ~ 0.7 (curve 3). 



The level of squeezing of the two modes is periodically 
repeated. The behavior of the two-mode variance also 
significantly depends on the phase matching condition, 
so that we may tune the phase sum to maximize the 
entanglement. We further choose for illustrations the 
phase sum corresponding to cos {T,9) — 1 for both opera- 
tional regimes. It seems perhaps more intuitive to study 
the time-behavior scaled on the parameter e — kE/j^, 
which contains the pump field amplitude and is an ad- 
justable parameter of the system. The dependence of V 
versus the scaled time-interval is shown in Fig. |21 where 
the three curves correspond to three different choices of 
the ratio e/x- Common to all curves is that the vari- 
ance is nonclassical and squeezed at least at the points 
of its minima imin, which can be obtained by the for- 
mula ctg (2^iniin) = s/ fJ-- In all cases, the maximal de- 
gree of two-mode squeezing Vmin = V (t) \t=t^i^ = 0.5 is 
achieved in the limit e —>■ x- For the curves on the Fig.|21 
the time-intervals corresponding to the minimal values 
of I4iin equal to eimin — 0.14 -I- 0.207rfc(curve 1), etmin — 
0.25 + 0.447rfc(curve 2), eimin ^ 0.39 + 0.987rA:(curve 3), 
(A: = 0,1, 2,...). 

If the opposite inequality holds, £ > Xj then the non- 
linear parametric interaction becomes dominant over the 
linear coupling and the variance is given by the following 
formula 



V{t) = 1 + ^ [cosh (2r;f ) - 1] 
V 

- - sinh (2ryi) cos (i;6') , 
V 



(60) 



X^- 



where 77 = ^ye^ 

For X ^ and cos (26*) — 1, from this formula we 
arrive at a well-known result, V {t) — exp(— 2ei), for 
two-mode squeezed state. This shows that in the limit of 
infinite squeezing the corresponding state approaches to a 




FIG. 3; Unitary time evolution of the variance V (t) for 
the range of e > x ^nd provided that cos (T,9) — 1. The 
parameters are: e/x ~ 1-1 (curve 1), e/x ~ 2 (curve 2) and 
e/x = 3 (curve 3) . 



simultaneous eigenstate of Xi — X2 and Yi +I2 , and thus 
becomes equivalent to the EPR state. Fig. O shows the 
behavior of V when the system operates in the regime e > 
X- This figure clearly shows that as the interaction time 
increases, the variance decreases and reaches its mini- 
mum. Then the squeeze variance exponentially increases 
with the growth of the interaction time. For the data 
in Fig. 131 we obtain that the time intervals for which the 
variance reaches the minima are etmin — 0.53(curve 1), 
einiin — 0.76 (curve 2), eimin — 0.93(curve 3). It is easy 
to check that these points of minima can be found by the 
formula ctgh{2r]tynin) — e/rj, provided that cos (£6*) — 1. 
We also conclude that the variance squeezes up to a 
certain interaction time, if e > %. It should be noted that 
although the time evolution of the variances are quite 
different for each of the operational regimes, the minimal 
values of the variance are described by the formula which 
is the same for both regimes 



14>i„ = 1^(t)k=t„„ = 



X 



X 



(61) 



With increasing e/x, in the regime e > Xi the minimal 
value of the variance decreases as Vmin ^ x/s ^ 1: which 
means that perfect squeezing takes place in the limit of 
infinite pump field. This result is not at all trivial for 
the system considered, even in the absence of dissipation 
and cavity-induced feedback, because the insertion of po- 
larization mixer usually destroys the two-mode squeezing 
produced by nondegenerate parametric down-conversion. 
The reason is that two-mode squeezed vacuum state is a 
supperposition of two-photon Fock states \n)-^ |n)2 and 
the polarization mixer destroys the Fock states having 
the same number of photons, i.e., ni = n2 = n. The de- 
tailed analysis of this problem can be found, for example, 
in [l|. 

Having discussed the CV entanglement for non- 
dissipative dynamics we now turn our attention to the 




FIG. 4: The minimized variance V versus the dimension- 
less amplitude of the pump field E/Eth ~ sjsth for both 
operational regimes. Sth and Em are given by the Eas. l25L 
12911 . The parameters are: xH = 0.1 (weak coupling), A/7 = 
10(curve 1); xH ~ 0.5(moderate coupling), A/7 = 3(curve 
2) and xll ~ 0.5, A/7 = 1 (curve 3). 



generation of entangled states of light beams in self-phase 
locked NOPO completely taking into the consideration 
dissipation and cavity-induced feedback. The results will 
be presented in a steady-state regime of generation. 



B. Entanglement in the self-phase locked NOPO: 
Sub-threshold regime 

Using formulas (|38|l . H55|l . after some algebra, we ob- 
tain the minimal variance for the steady-state regime in 
the following form 



V = 1 



eS^ 



725-4 -H A2 (5*2 - 2x2)^ 



5-4 _ 4A2;^2 



(62) 



It is easy to check that in the limit x — ^ the vari- 
ance coincides with the analogous one for the ordinary 
NOPO, V ^ V- = V^ ^ \ - e/{e -^ ^7 +^^)- We 
see that the minimal variance remains less than unity for 
all values of pump intensity and is a monotonically de- 
creasing function of e^ . For all parameters the maximal 
degree of two- mode squeezing y ~ 0, 5 is achieved within 
the threshold range. It is also easy to check that this 
expression is well-defined for all values of e^, including 
the vicinity of the threshold. One should keep in mind, 
however, that the linear approach used does not describe 
the threshold range where the level of quantum fluctua- 
tions increases. As a consequence, some matrix elements 
of (|38|l . (|48|) . H49|) increase infinitely in the vicinity of 
threshold. Nevertheless, it follows from 1)621) and further 
results H64|I - (I67|I . that such infinite terms are cancelled in 
the variance V , as well as in V_ , V+ for both operational 
regimes. This is not surprising, since such cancellation 
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of infinities in the quadrature amplitude variances takes 
place for the ordinary NOPO also. 

One of the differences between the squeezing effects 
of the ordinary and self-phase locked NOPO is that the 
variances V- and V+ for the ordinary NOPO are equal 



to each other, while for the self-phase locked NOPO they 
are in general different. The values of the non symmetric 
variances are expressed through R according to formula 
(|51|l . This quantity can be calculated with the help of 
the formula itsHl. The result is found to be 



J 



R 



exA 



5-4 _ 4A2^2 



725*4 + A2 (52 - 2x^f 

I 



2e 



(63) 



As we see, the final expressions (jHSJi l|63ll below the 
threshold are rather unwieldy. We show the correspond- 
ing numerical results on the Figs. ^|S1|H1 for illustration. 



C. Entanglement in the self-phase locked NOPO: 
Above-threshold regime 

Performing calculations for each of the cases A > 
and A < 0, we find the variance V in the above-threshold 
regime in the following form 

3 

4 



V 



1 



X 



l^l + lWTf ((£M.)'-l) 



41 A 



(64) 



To rewrite this expression through the original param- 
eters we should take into account that E/Eth = s/sth 
and £4^/7^ = (x/7 - |A| /jf + 1. The results ^ 
and H64|l for both operational regimes are summarized 
in Fig. 21 where the variance V is plotted as a function 
of the amplitude of the pump field. One can immedi- 
ately grasp from the figure that the sum of variances 

V = ^ {V+ + V-) remains less than unity for all nonzero 
values of the pump field, provided that x < |A|. This 
shows the nonseparability of the generated state. The 
maximal degree of entanglement is achieved in the vicin- 
ity of threshold, T^ ~ 0,5, if x/ |A| <C 1. In far above 
the threshold, E » Eth, V increases with mean photon 
numbers of the modes and reaches the asymptotic value 

V = 3/4 + x/4 |A|. It should also be mentioned that the 
result l|M|) is expressed through the scaled pump field 
amplitude e — kE/j^ and hence depends on coupling 
constants k and x- 

We stress, that the level of two-mode squeezing in the 
proposed scheme is limited due to the dissipation and 
pump depletion effects and reaches only 50% relative to 
the level of vacuum fluctuations if the pump intensity 
draw near to the generation threshold. Such limitation 
on the level of intracavity two- mode squeezing takes place 
also for the ordinary NOPO as was been demonstrated 
in Refs. [13, 1231 by analysis of quantum fluctuations in 
the near-threshold range. Analogous conclusion has been 
made for one- mode squeezed light generated in OPO |23l |. 



However, remind that these results pertain to the full 
squeezing of intracavity variances and not the spectra of 
the squeezing of output flelds. 

As the last step in connecting two-mode squeezing to 
observables of CV entanglement, we report the expres- 
sions of the non symmetric variances calculated with the 
help of the formulas (|5T|l . (|55l ) and (jSHJ. Upon evaluat- 
ing all required expectation values, we obtain 



R^ 



sgn{A) \A\-x 



1 



A 



V^TW^^Uh' 



(65) 



and hence 



V+ 



3 \ A \ +x + sgn (A) cos {A0) (| A | -%) 
1 + sgn (A) cos (A^) 



A^TTW^^Uh' 



(66) 



V- 



3\ A\+x- sgn (A) cos {Ad)i\A\~x) 

4|A| 
1 - sgn (A) cos {A9) 



4^l + {s^-ei,)h^' 



(67) 



For the case of A6 = 0, when the variances are maximally 
different, the results are reduced to 



V- = 



1 



X 

2A' 



V+ = 1- 



2^1 + {stk/lf {is/et,y - 1) 



(68) 



for A < 0. The case of A > is obtained from (|68|l by 
exchanging V+ — * V- and V- -^ V+. These results for 
both operational regimes are depicted on Fig. [S] Some 
features are immediately evident. First of all one can 
see that both variances are minimal in the critical range 
but show quite diflFerent dependences on the ratio e/eth 
above the threshold. Attentive reader may ask about de- 
pendence of the results obtained on parametric coupling 
constant k. We note in this connection that the squeezed 



11 




v.v 



FIG. 5: The variances V+ (curve 1) and V- (curve 2) versus 
the dimensionless amphtude of the pump field E/Eth = s/sth 
for the parameters: x/7 = 0.5, A/7 = 3. 




FIG. 6: The product of the variances V+ and V- versus the 
dimensionless amplitude of the pump field E/Eth = £/sth for 
the parameters x/7 ~ 0.1, A/7 = 10 (curve 1) and x/l = 
0.5, A/7 = 1 (curve 2). 



variances are expressed through the scaled pump field 
amplitude e = kE/js and hence depend in general on 
both coupling constants k and x- 

In terms of demonstrating the CV strong EPR entan- 
glement, one has to apply another criterion V+K_ < j. 
For the case of the symmetric uncertainties (V+ = V- = 
V), the product of the variances V"^ > j and hence the 
strong EPR entanglement can not be realized. In the 
general case, the product of the variances reads as 



V+V- = 1/2 _ ^2 ^Qg2 (-^gi) 



(69) 



It seems that V+V- lies belovif 1/4 at least for the relative 
phase Ad — ±TTm, (m — 1,2, ...), and in the vicinity of 
threshold, where F ~ 0, 5. However, for such selection of 
the phases we arrive at 



V+V- = 



A 



/ \ 



41 A 



^1 + [stHhf {{e/sthf - 

(70) 
It is easy to check that the product of the variances ex- 
ceeds J even in the vicinity of the threshold. This quan- 
tity for both operational regimes is illustrated in Fig. 
It should be noted again that a detailed analysis of this 
problem, must include more accurate consideration of 
quantum fluctuations in the critical ranges. 



V. CONCLUSION 

Our work demonstrates the possibility of creation of 
CV entangled states of light beams in the presence of 
phase localizing processes. We have called such quantum 
states the entangled self-phase locked states of light |21| 
and we have shown that they can be generated in NOPO 
recently realized in the experiment G*. This device is 



based on the type-II phase-matched down-conversion and 
additional phase localizing mechanisms stipulated by the 
intracavity waveplate. The novelty is that this device 
provides high level of phase coherence between the sub- 
harmonics in contrast to what happens in the case of ordi- 
nary NOPO, where the phase diffussion takes place. We 
have shown that both two mode squeezing and quantum 
phase locking phenomena are combined in such NOPO. 
This development paves the way towards the generation 
of bright CV entangled light beams with well-localized 
phases. It looks like that this scheme involving phase 
locking may be potentially useful for precise interferomet- 
ric measurements and quantum communications, because 
it combines quantum entanglement and stability of type 
- II phase matching with effective suppression of phase 
noise. Particularly, in practice, this phase locked scheme 
can effectively be used in homodyne detection setups, as 
well as for observation of an interference with nonclassi- 
cal light beams. We believe, that such source of bright 
entangled light providing phase coherence will also be 
applicable for realizing CV quantum teleportation. The 
price one has to pay for these advantages is the small ag- 
gravation of the degree of CV entanglement in self-phase 
locked NOPO in comparison with the case of ordinary 
NOPO. The quantum theory of self-phase locked NOPO 
has been developed in linear treatment of quantum fluc- 
tuations for both below- and above-threshold regimes of 
generation. We have studied the CV entanglement as 
two-mode squeezing and have shown that entanglement 
is present in the entire range of pump intensities. In all 
cases the maximal degree of two-mode squeezing V ~ 0, 5 
is achieved in the vicinity of the threshold, provided that 
the coupling constant x is much less than the mode de- 
tunings. We have demonstrated that, as a rule, highest 
degree of CV entanglement occurs for weak linear cou- 
pling strength x- It has also been shown that the amount 
of the entanglement can be controlled via the phase differ- 
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ence $^. The other peculiarities of the system of interest 
have been estabhshed for the case of unitary dynamics 
which may be reahzed for the short interaction times in 
the case of pulsed pump fields. One of these concerns the 
presence of two operational regimes generating two-mode 
squeezing. If the linear coupling between subharmonics 
dominates over the parametric down-conversion, e < Xj 
we have observed a periodic evolution of the squeezed 
variance. Maximal degree of squeezing has been 1/ ~ 0, 5 
in this regime. If the parametric interaction becomes 
dominant, £ > Xj the more high degree of two-mode 
squeezing can be obtained, Vmin = j^ < 0.5, but up 
to certain interaction time. 

In our analysis we have not investigated all possible 
quantum effects of self-phase locked parametric dynam- 
ics. In particular, we have noted that the system consid- 
ered displays different types of quantum correlations, but 



we have not analyzed their connection with all possible 
kinds of the entanglement. This analysis might involve 
also the case of unitary dynamics, where the generalized 
two-mode squeezed state has been presented. Consider- 
ation of quantum fiuctuations in the near threshold op- 
erational range of self-phase locked NOPO also deserves 
special attention for more accurate identification of CV 
entanglement. These topics are currently being explored 
and will be the subject of forthcoming work. 
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